Rayleigh's Stretched String 



Mark B. Villarino 
Depto. de Matematica, Universidad de Costa Rica, 
2060 San Jose, Costa Rica 

January 2, 2008 



Abstract 

We obtain a rigorous a priori upper and lower bounds to the exact period of the 
celebrated Rayleigh stretched-string differential equation. 

1 Introduction 

In his monumental treatise The Theory of Sound (see pages 45-46 [2]) Lord Rayleigh 
proposes the following problem: 

A horizontal elastic steel wire having length 2Lq and spring constant a is stretched 
to length 2L. A body, having a mass m much greater than the mass of the wire, 
is tied to the center of the wire and is put in motion on a vertical line through 
the center of the wire. The only forces on the body are produced by tension in the 
wire (there is no gravity, no damping). It is required to study the motion of the 
body. 

(We have used Agnew's formulation [1], p. 17-18.) Let y{t) = yhe the vertical displacement 
of the mass m in time t, where y is positive, 0, or negative according as the body is above, 
or at, or below the wire in equilibrium or neutral position. 

By the law of HoOKE, the magnitude of the elastic force, or the tension, T{y), exerted 
by each half of the stretched string is equal to: 
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and therefore the vertical component of the total force, exerted by both halves, which is the 
only component acting to move the mass m is given by 

vertical component of T(y) = —2a- ^ ~ ~ — ^ ■> (1-1) 

Lq a/ L2 + y2 



where we have employed our sign conventions. 

But, by Newton's second law of motion, that same force, with opposite sign, is given 
bym||. 

Therefore, equating the two forces, and bringing the expression for the elastic force on 
the same side as the force given by Newton's second law, we find that the algebraic sum 
of the forces is zero and dividing by m, we conclude that the equation of motion of the mass 
m is given by: 
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This differential equation (II. 2p is nonlinear and very complicated, and we have never seen 
an exact treatment of it in the literature. 

Rayleigh, himself, simplifies the equation by the following physical reasoning; 

"The tension of the string in the position of equilibrium depends on the 
amount of the stretching to which it has been subjected. In any other position 
the tension is greater; but we limit ourselves to the case of vibrations so small that 
the additional stretching is a negligible fraction of the whole. On this condition, 
the tension may be treated as a constant." 

(Our italics. )We just saw that the magnitude of the tension, T{y), (or the force) from 
each half of the string on m is given by 
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provided the stretching is well within elastic limits. So, Rayleigh effectively takes y = in 
(11. 3p so that the equation (11. 2p becomes 



d'^y 2a L — Lq 
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or, using the notation (O]) (with y = 0, i.e., T := T(0) = cr^^), 

2T 

y + —y = 
mL 
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where the "dots" refer to time derivatives. If we assume that the initial displacement at time 
t = is yo, then the solution to (11. 5p is 
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which is simple harmonic motion with period P given by: 
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Rayleigh does not discuss how close his approximate period P ( II. 7p is to the exact 
period P. Indeed, we have been unable to to find any error analysis of Rayleigh's solution 
in the literature. Yet, such a famous and classical differential equation merits an investigation 
into the accuracy of its approximate solution. 

Therefore we offer the following theorem to fill this gap. 

Theorem 1. Let P be the true period of oscilation of the mass m. Then Rayleigh's 
approximate period P overestimates the true period, P. Indeed, the following inequalities 
are valid: 

P _ 
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Our lower bound shows the dependence of the period on the reciprocal, of the initial 
displacement, and tends to the upper bound as Hq tends to 0. 

2 The Formula for the Period 



The differential equation (11. 2p does not have the variable t nor the first derivative, ij , ap- 
pearing explicity. Therefore, we can transform it into a separable equation for the velocity, 
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Substituting the right-hand side of (12. 2p into (II. 2p . transposing the expression for the tension 
to the right-hand side, separating the variables and forming the indefinite integral of both 
sides we obtain 



(2.3) 
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where C is the constant of integration. At time t = 0, the conditions are y{0) = yo and 
v{0) = 0. Therefore, evaluating the constant C, recalling that w = taking the square 
root of both sides, and using a negative sign because the displacement, y{t), is a decreasing 
function during the first quarter oscilation, we obtain: 
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Multiplying the quarter-period by 4, we obtain the following formula for the full period, P, 
of oscilation of the mass: 



3 



Theorem 2. The true period, P, of oscillation of the mass m is given by: 




This integral is very comphcated for exact computation. 
Indeed, if we make the change of variable 



we obtain the formula 
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P = P{zo)-P{L). 



(2.5) 
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The integral in (12. 9p is a very complicated algebraic combination of an elliptic integral of the 
third kind and an elliptic integral of the first kind, and we will not pursue this line of inquiry 
any further. 

Instead, we turn our attention to the upper and lower bounds for the true period, P. 



3 Upper and Lower Bounds for the True Period 

Proof of the bounds on the true period. We begin with the exact formula (13.11) for the true 
period, P: 
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Taking y = = in the sum -^/Z^~+~^ + ^JL? ^y^ only makes the integrand larger 
and we conclude that 



P <4 
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and we have proved the validity of the upper hound: 




(3.2) 



For the lower bound, we proceed similarly, and we start by observing that if we take 
y = yo va. the sum ^ I? ^ y^ + \/ L'^ ^ y\ the integrand becomes smaller and we conclude 
that 
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where, twice, we have used the inequality 
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Therefore, we have proved the validity of the lower bound: 




This completes the proof of the theorem 11.81 
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Corollary 1. If P denotes Rayleigh's approximate period {\1.7\i . then, the relative error 

P-P 



R :-- 
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in the approximation P ^ P, satisfies the inequality 
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This shows that the lower bound to the relative error is proportional to the square of the 
initial displacement and tends to zero with a "quadratic" rate of convergence. 
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